ADP-98-48/T317, OSUNT98-13 



Charmed Mesic Nuclei: 
Bound D and D states with 208 Pb 

K. Tsushima 1 , D.H. Lu 1 , A.W. Thomas 1 , K. Saito 2 and R.H. Landau 3 
1 Department of Physics and Mathematical Physics, and Special Research Center for the 



(N 
> 

o 
o 



Subatomic Structure of Matter, University of Adelaide, Adelaide, 5005 Australia 



2 Physics Division, Tohoku College of Pharmacy, Sendai 981-8558, Japan 



3 Department of Physics, Oregon State University, Corvallis, OR 97331, USA 

(February 9, 2008) 

We show that the D~ meson will form narrow bound states with 208 Pb. Mean field 
potentials for the D°, D° and D~ in 208 Pb are calculated self-consistently using the 

00 

quark-meson coupling (QMC) model in local density approximation. The meson- 
208 Pb bound state energies are then calculated by solving the Klein-Gordon equation 
^ ' with these potentials. The experimental confirmation and comparison with the D° 

a' 

and D° will provide distinctive information on the nature of the interaction between 



the charmed meson and matter. 
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PACS numbers: 21.30.Fe, 24.10.Jv, 14.40.Lb, 21.10.Dr, 11.30.Rd 



In relativistic models of nuclear structure like QHD or QMC p-4T0f, the isoscalar-scalar 
meson, a, is responsible for a large reduction in the mass of the nucleon (m^ — » m* N ). Although 
there is no firm evidence of a a meson with mass in the range 500 - 600 MeV (of. the Particle Data 
Group discussion of the /o(400 — 1200) and a |ll]]), there is considerable empirical justification 
for using the a as a phenomenological representation of correlated two-pion-exchange between 



nucleons [12]. In QMC the justification for coupling the a to the confined light quarks (hereafter 
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referred to as q) is that dynamical chiral symmetry breaking requires that the confined quarks 
couple to pions. Quarks in different hadrons can interact by exchanging two pions and this is 
represented phenomenologically by cr-exchange. Using this model it is possible to investigate the 
reduction of the masses not only of baryons, but also of mesons containing light quarks (See 



Refs. [13-jn| for other approaches not all of which result in a mass reduction for the mesons.) 

The result for the w-meson within relativistic mean field models is especially interesting. 
Because it consists of almost pure q-q pairs (ideal mixing) , in an isoscalar nucleus the q and q feel 
equal and opposite vector potentials. This means that the effect of the mean cr-field is unmasked 
which leads to expectations of quite deeply bound w-nucleus states [f7,14,20, ^l|]23}l . These states 



are currently the subject of intense experimental investigation, with the most promising involving 
recoilless production in the (d, 3 He) reaction at GSI [20,21|. 



In this paper we consider a possibility that is in some ways even more exciting, in that 
it promises more specific information on the relativistic mean fields in nuclei and the nature of 
dynamical chiral symmetry breaking. We focus on systems containing an anti-charm quark and a 
light quark (cq), which have no strong decay channels if bound. If we assume that dynamical chiral 
symmetry breaking is the same for the light quark in the charmed meson as in purely light-quark 
systems, we expect the same a-q coupling constant. (On the one hand, this assumption seems 
quite reasonable, given the success in reproducing the pion coupling constants for various charmed 



baryons using a model based on a common light quark pion coupling [22], while on the other hand, 
the present predictions could also be viewed as an independent test of this picture.) In the absence 
of any strong interaction, the D~ will form atomic states, bound by the Coulomb potential. We 
use the QMC model to estimate the effect of the strong interaction. The resulting binding for, say, 
the Is level in 208 Pb is between ten and thirty MeV and should provide a very clear experimental 
signature. 

Systems of the form qc are also extremely interesting because the coupling of the mean vector 
field to the light anti-quark will be attractive. Indeed, we expect a D meson (cq) to experience an 
attraction in excess of 100 MeV in an average size atomic nucleus. Unfortunately, the D meson 
in matter will also couple strongly to open channels such as DN — ► B c (tt's), with B c a charmed 



baryon. Our present knowledge does not permit an accurate calculation of the corresponding 
widths, which may be 10-100 MeV. At the lower end of this range, such states should be able to 
be seen as very deeply bound mesic-nuclei, while at the upper end they may not be detectable. 
We regard the widths as an experimental issue at present and show the Pb bound state energies 
without the effect of absorption. 

A detailed description of the Lagrangian density and the mean-field equations of motion 
needed to describe a finite nucleus is given in Refs. At position r in a nucleus (the coordinate 

origin is taken at the center of the nucleus), the Dirac equations for the quarks and antiquarks in 
the D and D meson bags are given by 



0. 



V 



^d(x) 

[*7 ' 9 X ~ m c ] ip c (x) (or ip 5 (x)) = 0. 
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(1) 

(2) 
(3) 



The mean- field potentials for a bag centered at position f in the nucleus are defined by V§(r) = 
g%cr(r),VJ(r) = g^u>(f) and V!jf(r) = g%b(r), with g%,g% and g q p the corresponding quark and 
meson-field coupling constants. (Note that we have neglected a possible, very slight variation of 
the scalar and vector mean-fields inside the meson bag due to its finite size ||.) The mean meson 
fields are calculated self-consistently by solving Eqs. (23) - (30) of Ref. [||, namely, by solving a 
set of coupled non-linear differential equations for static, spherically symmetric nuclei, resulting 
from the variation of the effective Lagrangian density involving the quark degrees of freedom and 
the scalar, vector and Coulomb fields in mean field approximation. 

The normalized, static solution for the ground state quarks or antiquarks in the meson bags 
may be written as: 



$ f {x) = N f e- i£ ^ R H f (x), (j = D, D), 



(4) 



where f = u,u, d, d, c, c refers to quark flavors, and Nf and ipf(x) are the normalization factor and 
corresponding spin and spatial part of the wave function. The bag radius in medium, R*, which 
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depends on the hadron species to which the quarks and antiquarks belong, will be determined 
through the stability condition for the (in-medium) mass of the meson against the variation of the 
bag radius ||||j6],[j]] (see also Eq. @). The eigenenergies e/ in Eq. (Q) in units of 1/Rj are given by 



( e d (r) ^ 



n* q (r)±RUV«(r) + -V«(r 



n* q (r)±R*[V«(r)--V«(r)), 



J 



(5) 

(6) 
(7) 



where 0*(r) = ^Jx 2 + (R^m*) 2 , with m* = m q — g$cr(r) and f2 c (r) = Jx 2 + {R*m c ) 2 . The bag 
eigenfrequencies, x q and x c , are determined by the usual, linear boundary condition HQ. (Note 
that the lowest eigenenergy value for the Dirac equation (Hamiltonian) for the quark, which is 
positive, should be identified with a constituent quark mass.) 

The D and D meson masses in the nucleus at position r, are calculated by: 



n*(r) + n c (r) - z< 4 
m*{r) = ^ W K> 3 + -irRfB, 



dm*j (r) 



dRi 



0, 



(j = D,D). 



(8) 
(9) 



In Eq. (|8j), the Zj parametrize the sum of the center-of-mass and gluon fluctuation effects, and are 
assumed to be independent of density. The parameters are determined in free space to reproduce 
their physical masses. 

In this study we chose the values m q = m u = = 5 MeV and m c = 1300 MeV for the 
current quark masses, and Rn = 0.8 fm for the bag radius of the nucleon in free space. Other 
input parameters and some of the quantities calculated are listed in Table |. The parameters at 
the hadronic level associated with the core nucleus can be found in Refs. We stress that 

while the model has a number of parameters, only three of them, g%, and g^, are adjusted to 
fit nuclear data - namely the saturation energy and density of symmetric nuclear matter and the 
bulk symmetry energy. None of the results for nuclear properties depend strongly on the choice 
of the other parameters - for example, the relatively weak dependence of the final results for the 



properties of finite nuclei, on the chosen values of the current quark mass and bag radius, is shown 
explicitly in Refs. Exactly the same coupling constants, g%, and g^, are used for the light 

quarks in the mesons as in the nucleon. However, in studies of the kaon system, we found that it 
was phenomenologically necessary to increase the strength of the vector coupling to the non-strange 
quarks in the K + (by a factor of 1.4?) in order to reproduce the empirically extracted if + -nucleus 
interaction ||. It is not yet clear whether this is a specific property of the K + , which is a pseudo- 
Goldstone boson, or a general feature of the interaction of a light quark, perhaps associated with 
the Pauli exclusion principle. In any case, we show results for the D binding energies with both 
choices for this potential, in order to test the theoretical uncertainty. The u mean field potential 
with the larger coupling will be labelled V£ (= 1A 2 V^). 

Through Eqs. (||) - @ we self-consistently calculate effective masses, m*(f) (j = D,D), and 
mean field potentials, V§ w „{f), at position r in the nucleus. The scalar and vector potentials felt 
by the hadron, which will depend only on the distance from the center of the nucleus, r = |r|, are 
given by: 

V s j (r) =m*(r) -nij, (10) 
Vr(r) = V2(r)-±V*(T)-A(r), (11) 
vf\r) = VZ{r)+ l -V«(r), (12) 
K D °W = -(^(r) + iy;(r)), (13) 

where A(r) is the Coulomb interaction between the meson and the nucleus. Note that the p meson 
mean field potential, VS(r), is negative in a nucleus with a neutron excess, such as e.g., 208 Pb. For 
the larger u> meson coupling, suggested by K + A scattering, VJ(r) is replaced by VJ(r). 

Before showing the calculated potentials for the D~ in 208 Pb which is particularly interesting 
in view of the strong Coulomb field, we first show in Fig. [I] the mass shift of the D(D) meson, 
calculated in symmetric nuclear matter (no contributions from the p and Coulomb fields). As is 
expected, the relation (m D jj — m* D £,) ~ |(mjv — m* N ) is well realized p. 

Next, in Fig. || we show the sum of the potentials for the D in 208 Pb for the two choices of 
V S D (r) + V® (r) (the dashed line corresponds to VJ(r) and the dotted line to VJ(r)). Because 
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the D~ meson is heavy and may be described well in the (nonrelativistic) Schrodinger equation, 
one expects the existence of the ^Pb bound states just from inspection of the naive sum of the 
potentials, in a way which does not distinguish the Lorentz vector or scalar character. 

Now we are in a position to calculate the bound state energies for the D and D in nuclei, 
using the potentials calculated in QMC. There are several variants of the dynamical equation for 
a bound meson- nucleus system. Consistent with the mean field picture of QMC, we actually solve 
the Klein-Gordon equation: 

[V 2 + (E j - VJ(r)) 2 ~ mf (r)] ^(r) = (14) 

where Ej is the total energy of the meson (the binding energy is Ej — raj). To deal with the long 
range Coulomb potential, we first expand the quadratic term (the zeroth component of Lorentz 
vector) as , (Ej-Vj(r)) 2 = E] + A 2 (r) + V 2 p (r) + 2A(r)V ujp (r) - 2E J (A(r) + V up (r)), where V up (r) 
is the combined potential due to oj and p mesons {V up (r) = V£(r) — \V p q {r) for D~). Then Eq. fli~4| ) 
can be rewritten as an effective Schrodinger-like equation, 

V 2 



2rrij 



+ Vj(Ej,r) 



$i(r) = ''^n!/' 1 ^^ (15) 



where <&j(r) = 2nij(j)j(r) and Vj(Ej, r) is an effective energy-dependent potential which can be split 
into three pieces (Coulomb, vector and scalar parts), 

JV J ' 1 rrij y ' 2mj 2mj v ' 

Note that only the first term in this equation is a long range interaction and thus needs special 



treatment, the second and third terms are short range interactions. In practice, Equation (15) is 
first converted into momentum space representation via a Fourier transformation and is then solved 
using the Kwon-Tabakin-Lande technique |24|] . We would like to emphasize that no reduction has 
been made to derive the Schrodinger-like equation, so that all relativistic corrections are included 
in our calculation. The calculated meson-nucleus bound state energies for 208 Pb, are listed in 
Table 0. 

The results show that both the D~ and D° are bound in 208 Pb with the usual oj coupling 
constant. For the D~ the Coulomb force provides roughly 24 MeV of binding for the Is state, and 
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is strong enough to bind the system even with the much more repulsive uj coupling (viz., 1A 2 V£). 
The D° with the stronger uj coupling is not bound. Note that the difference between D° and D~ 
without the Coulomb force is due to the interaction with the p meson, which is attractive for the 
D° but repulsive for the D~. For completeness, we also calculated the binding energies for the D°, 
which is deeply bound since the u interaction with the light antiquarks is attractive. However, the 
expected large width associated with strong absorption may render it experimentally inaccessible. 
It is an extremely important experimental challenge to see whether it can be detected. 

The eigenfunctions for the Schrodinger-like equation are shown in Fig. ||, together with the 
baryon density distribution in 208 Pb. For the usual ui coupling, the eigenstates (Is and lp) are well 
within the nucleus, and behave as expected at the origin. For the stronger u coupling, however, 
the D~ meson is considerably pushed out of the nucleus. In this case, the bound state (an atomic 
state) is formed solely due to the Coulomb force. An experimental determination of whether this 
is a nuclear state or an atomic state would give a strong constraint on the u coupling. We note, 
however, that because it is very difficult to produce .D-mesic nuclei with small momentum transfer, 
and the D-meson production cross section is small compared with the background from other 



channels, it will be a challenging task to detect such bound states experimentally [25|. 

To summarize, we have calculated the Pb, Pb and |^Pb bound state energies in QMC. 
The potentials for the mesons were calculated self-consistently in local density approximation. In 
spite of possible model-dependent uncertainties, our results suggest that the D~ meson should be 
bound in 208 Pb due to two quite different mechanisms, namely, the scalar and attractive a mean 
field even without the assistance of the Coulomb force in the case of the normal vector potential 
(V£), and solely due to the Coulomb force in the case of the stronger vector potential (VJ). (We 
recall that the kaon is a pseudo-Goldstone boson and expected to be difficult to treat properly with 
the usual bag model. Thus, the analysis of Ref. || on the vector potential for the light quarks 
inside the kaon bag may not be applicable to the light quarks inside the D-meson.) 

Thus, whether or not the L> - 208 Pb bound states exist would give new information as to 
whether the interactions of light quarks in a heavy meson are the same as those in a nucleoli. The 
enormous difference between the binding energies of the D° (~ 100 MeV) and the D° (~ 10 MeV) 



is a simple consequence of the presence of a strong Lorentz vector mean-field, while the existence of 
any binding at all would give us important information concerning the role of the Lorentz scalar a 
field (and hence dynamical symmetry breaking) in heavy quark systems. In spite of the perceived 
experimental difficulties, we feel that the search for these bound systems should have a very high 
priority. 
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TABLES 

TABLE I. The physical masses fitted in free space, the bag parameters, z, and the bag radii in 
free space, R. The quantities with an asterisk, *, are those quantities calculated at normal nuclear 
matter density, po = 0.15 fm™ 3 . They are obtained with the bag constant, B = (170.0 MeV) 4 , 
current quark masses, m u = rrid = 5 MeV and m c = 1300 MeV. (Note that the free space widths 
for D and D mesons are negligible pT| .) 





mass (MeV) 


z 


R (fm) 


m* (MeV) 


R* (fm) 


N 


939.0 (input) 


3.295 


0.800 (input) 


754.5 


0.786 


D,D 


1866.9 (input) 


1.389 


0.731 


1804.9 


0.730 



TABLE II. Calculated D~ , D° and D° meson bound state energies (in MeV) in 208 Pb for 
different potentials. The widths for the mesons are all set to zero, both in free space and inside 
208 Pb. Note that the D° bound states energies calculated with VJ will be much larger than those 



calculated with VJ (in absolute value). 



state 




d-{V3) 


D-(Vj, no Coulomb) 


&>{V$) 




D°(VJ) 


Is 


-10.6 


-35.2 


-11.2 


unbound 


-25.4 


-96.2 


IP 


-10.2 


-32.1 


-10.0 


unbound 


-23.1 


-93.0 


2s 


-7.7 


-30.0 


-6.6 


unbound 


-19.7 


-88.5 
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FIGURES 



t — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

FIG. 1. Shift of the mass of the nucleoli and the D or D meson. (Normal nuclear matter 
density, po, is 0.15 fm~ 3 .) 
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FIG. 2. Sum of the scalar, vector and Coulomb potentials for the D~ meson in 
208 Pb for two cases, (m* D _(r) - m D -) + V£(r) + \V!}{r) — A(r) (the dashed line) and 
(m* D _(r) - m D -) + FJ(r) + \Vj}(r) - A{r) (the dotted line), where FJ(r) = 1.4 2 FJ(r). 
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FIG. 3. The Schrodinger-like bound state wave functions of the D~ meson in 208 Pb, for two 
different to meson coupling strengths. See also the caption of Fig. |2[ The wavefunction is normalized 
as follows: J °°(ir 47rr 2 |$(r)| 2 = 1. 
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